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ABSTRACT 

Local Group dwarf spheroidal satellite galaxies are the faintest extragalactic stellar systems known. 
We examine recent data for these objects in the plane of the Baryonic Tully-Fisher Relation (BTFR). 
While some dwarf spheroidals adhere to the BTFR, others deviate substantially. We examine the 
residuals from the BTFR and find that they are not random. The residuals correlate with luminosity, 
size, metallicity, ellipticity, and susceptibility of the dwarfs to tidal disruption in the sense that 
fainter, more elliptical, and tidally more susceptible dwarfs deviate farther from the BTFR. These 
correlations disfavor stochastic processes and suggest a role for tidal effects. We identify a test to 
distinguish between ACDM and MOND based on the orbits of the dwarf satellites of the Milky Way 
and how stars are lost from them. 

Subject headings: dark matter — galaxies: dwarf — galaxies: formation — galaxies: halos — Local 
Group 



1. INTRODUCTION 

Recent years have seen enormous progress in the 
discovery and measurement of the tiny dwarf satel- 
lite galaxies of the Local Group. These include the 
long known, 



"classical" dwarf spheroidal galaxies ([Mate* 
1998) as well as the more recently dis covered "ultra- 



faint" dwarfs and the satellites of M31 (Willma n et al 
20051: iZucker et al.ll200d:lGrillmairll2006tlMaiewski et al 
2007tlBelokurov et al.ll2007HGrillmairl2009HMartin et al 
2009). In addition to identifying these systems, kine- 



matic data from measuring the velocities of individual 
stars has become available for many systems. These now 
consist of t housands of indivi dual stars for the classi- 
cal dwarfs (| Walker et al.l 120071 ) . wi th rapidly improving 
data for the other t ypes of systems (|Sirnon fc Gehall2007l : 
IKalirai et ani2010h . 

The Local Group dwarfs appear to be the most 
dark matter dominated objects in the universe 
dMatedll998HWilkinson et al.ll20q2t iGilmore et al.1120071: 
Simon fc Gehal 120071: iKqch et al.l l2007bt iStrigari et al.1 
20081 iWalker et al.l l2009t iWolf et al.l 120101 ) . consistent 
with the trend of dark matter do mination increasing 
with decreasing surface brig htness (|de Blok et al.|[T996l : 
iMcGaugh fc de Blold [l998af As such, they provide a 
unique probe of structure formation at the smallest ac- 
cessible scales. They presumably reside in the sub-halos 
thought to inhabit the large dark matter halos of galax- 
ies like the Milky Way (MW) and M31. If so, physical 
processes specific to their environment, like tidal disrup- 
tion and ram pressure stripping, may play a role in the 
evolution of the luminous content of the dwarfs as they 
orbit the primary structure. It is therefore interesting to 
investigate whether, and the extent to which, the dwarfs 
obey the scaling relations established for brighter galax- 
ies. 

In this paper we investigate how the dwarfs behave 



in the Baryonic Tully-Fisher plane. Rotating disk 
galaxies define a tight relation between baryonic mass 
and rotation velocity (jMcGaugh et al.l 120001: iVerheiien 



2001). This has recently bee n extended (jStark et al 



2009; iTrachternach et aTl [2009) to low rotation veloci- 
ties (~ 20 kms -1 ) comparable to the Local Group dwarf 
spheroidals. An obvious question is whether the Local 
Group dwarfs that are satellites of the MW and M31 
continue this relation. 

2. DATA 



We adopt the data compilation of IWolf et al.l (|2010l ) 
for the MW satellites au gmented by the rec ent data for 
M31 dwarfs presented bv lKalirai et al.l (I2010D. These are 
small , pressure supported systems (|Mateolll99 8: Strigari 
l2010f ) whose baryonic mass is in the form of stars Q In 
order to place them on the same Baryonic Tully-Fisher 
Relation (BTFR) as rotating disks we need to estimate 
stellar mass and rotation velocity in a consistent manner. 

The stellar populations of the Local Group dwarfs are 
not simply ancient single burst populations. While some 
are old, others show clear signs of continuing intermit- 
tent s t ar forming epi s odes j Grcbcl 1997; Hernande z et al.l 
I2000D . iMateo et all (|1998fl have estimated the mass-to- 
light ratios Tjf of the stars of the classical d warfs from re- 
solved observations of their stellar content. iMartin et al.l 
(|2008aj) have done the same for many of the ultrafaint 
dwarfs. We adopt their stellar mass-to- light ratios for 
the objects studied. The mean of each of these samples 
is 1.3 Mf?> /Lfn if we adopt the Kroupa IMF in the tab- 
ulation of lMartinet al. (2008a). For objects not specifi- 
cally addressed by these studies, namely Leo T and the 
satellites of M31, we use the mean mass-to-light ratio. 

Leo T is the one exception to this, with a non-negligible amount 
of HI gas (Grccvich & Putman 2009). We include the gas mass in 
the baryonic sum. 
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Table 1 

Baryonic Tully-Fisher Data and Residuals 



Dwarf 



D 



Vc 



T 



v 



log M b 



logF 6 [Fe/H] <5[Fc/H] Refs. 



1 

11 



Carina 
Draco 
Fornax 
Leo I 
Leo II 
Sculptor 
Sextans 
Ursa Minor 
Bootes I 
Canes Ven. 
Canes Ven. 
Coma Ber. 
Hercules 
Leo IV 
Leo T 
Segue 1 
Ursa Major I 
Ursa Major II 
Willman 1 
And I 
And II 
And III 
And VII 
And X 
And XIV 



101 ± 5 
76 ± 5 
138 ± 8 
250 ± 30 
205 ± 12 
79 ± 4 
86 ± 4 
66 ± 3 
66 ± 3 
218 ± 10 
160 ± 4 
45 ± 4 
132 ± 12 
160 ± 15 
407 ± 38 
28 ± 2 
97 ± 4 
36 ± 5 
43 ± 7 
57 ±24 
194 ± 18 
75 ± 24 
215 ±35 
109 ± 37 
186 ± 87 



11.1 ± 0.4 
17.5 ±0.9 
18.5 ±0.3 
15.5 ±0.7 

11.4 ±0.8 

15.5 ±0.3 

12.3 ±0.5 

20.0 ± 1.1 
15.7 ± 3.9 

13.2 ±0.8 
7.9 ± 1.8 

8.0 ± 1.4 
8.9 ± 1.6 

5.8 ±3.0 

13.5 ± 2.8 
7.4 ± 2.0 

13.1 ± 1.8 

11.7 ± 2.5 

6.9 ± 1.6 

18.4 ± 1.9 

12.6 ± 1.4 

8.1 ± 3.1 

16.8 ± 2.8 
6.8 ±2.1 
9.4 ± 2.3 



1.0 
1.5 
1.1 
0.9 
1.4 
1.4 
1.3 
1.5 
1.1 
1.3 
1.0 
1.3 
1.0 
1.0 
1.3 
1.8 
1.4 
1.4 
1.5 
1.3 
1.3 
1.3 
1.3 
1.3 
1.3 



64 ±0.41 
50 ± 0.42 
27 ± 0.42 

65 ±0.42 
03 ± 0.42 
56 ± 0.42 

89 ± 0.42 

76 ± 0.43 
49 ± 0.41 
48 ± 0.41 

90 ± 0.45 
68 ± 0.45 
05 ± 0.43 
94 ± 0.46 
92 ± 0.40 
78 ± 0.50 
29 ± 0.42 
75 ± 0.45 
18 ±0.50 

77 ± 0.40 
09 ± 0.41 
11 ± 0.42 
37 ± 0.42 
29 ± 0.45 
37 ± 0.45 



-0.18 
-1.12 
0.55 
0.24 
0.16 
0.14 
-0.13 
-1.09 
-1.94 
-0.66 
-1.35 
-1.58 
-1.39 
-0.75 
-0.25 
-2.35 
-1.83 
-2.17 
-1.84 
0.06 
1.03 
0.82 
0.82 
0.32 
-0.16 



-1.80 
-1.99 
-1.29 
-1.31 
-1.74 
-1.81 
-2.07 
-2.03 
-2.50 
-2.08 
-2.19 
-2.53 
-2.58 
-2.58 
-2.02 
-3.30 
-2.29 
-2.44 

-1.45 
-1.64 
-1.78 
-1.40 
-1.93 
-2.26 



0.30 
0.32 
0.46 
0.25 
0.23 
0.34 
0.36 
0.32 

0.46 
0.58 
0.45 
0.51 
0.75 
0.54 
0.50 
0.54 
0.57 

0.37 
0.34 
0.27 
0.30 
0.11 
0.31 



1,3,10,12 
1,4,10,12 
1,3,10,12 
1,5,10,12 
1,6,10,12 
1,3,10,12 
1,3,10,12 
1,4,10,12 
1,7,11,12 
1,8,11,12 
1,8,11,12 
1,8,11,12 
1,8,11,12 
1,8,11,12 
1,8,12 
1,9,11,12 
1,8,11,12 
1,8,11,12 
1,11,12 
2,12 
2,12 
2,12 
2,12 
2,12 
2,12 



Refe rences. — Distanc es and kinematic d ata: l.lWolf et al.| ||2 010D 2. |I^al^raiet {2 010I). Metallicity 
data: 3. [H"elmi et alJl | 2006P 4 . IWinnickl 112003) 5.IKoch et al.l !2007cl) ir iKoch et al.l d20p7al') 7.|Martin et all 
( 2007 ) 8. | Kirbvet al] \2003) 9. IGeha et al.l 120091) . Stellar mass-to-light ratios: 10. IMateo et all (jxQ98D 
ll. lMartin et all d2008al) . Baryonic mass and BTFR residuals: 12. this work. 

Note. — Columns: (1) Name of the dwarf satellite. (2) The distance of the dwarf from the center of its 
host in kpc. (3) The circular velocity (equation^ in km s — . (4) The V-band stellar mass-to-light ratio 
(in Mq/Lq) assumed to compute baryonic mass. (5) The logarithm of the baryonic mass (M^ — T^Ly) 
in Mq. (6) The logarithm of the deviation from the BTFR (equation [2} : F b = Af b /(45V C 4 ). (7) The 
mean metallicity and (8) the width of the metallicity distribution. (9) The source of the data. 



These stellar mass-to-light ratios are subject to consid- 
erable uncertainty, particularly in the IMF. We therefore 
adopt a conservative uncertainty in this quantity of 0.4 
dex. We also propagate the uncertainties in luminos- 
ity and distance, but in most cases the uncertainty in 
the mass is dominated by that in the mass-to-light ratio: 
M b = TYL v . 

The stellar mass is woefully inadequate to explain the 
observed velocity dispersions of these systems: they ap- 
pear to be dark matter dominated. We adopt the line of 
sight velocity dispersions oi os tabulated by IWolf et al 



(2010) for the Milky Way satellites and Kalirai et al 
(2010) for the faint satellites of M31. In order to re- 
late this value to the characteristic circular velocities of 
rotating disk galaxies, we assume 

V c = V3a lo8 . (1) 

We take this to be the circular velocity characteristic 
of the dark matter halo of each dwarf. This is not pre- 
cisely the same quantity as measur ed in disks. Geometry 
can ca use differences up to ~ 20% (|McGaugh fc de Blokl 
19983), m the sense that flattened systems of the same 
mass rotate faster. In addition, the radius at which oi os is 
measured is typically modest compared to that of spirals. 
Nevertheless, adopting the velocity dispersion measured 
at the 3D half light rad ius minimizes th e uncertainty 
due to orbital anisotropy ()Wolf et al.ll2010D and provides 
a lower limit on the characteristic circular velocity at 
larger radii. Moreover, the velocity dispersion profiles of 
the well-observed classical dwarfs tend to be rather flat 



(Walk er" et al.| [2007). so this is the best available estima- 
tor of the gravitational potential. The circular velocities 
calculated in this manner are listed in Table [TJ We also 
tabulate other useful quantities, computed here or col- 
lected from the literature, in Tables [1] and [5J 

3. INTERPRETATIONS 

Provided with estimates of the baryonic masses and 
circular velocities of the Local Group dwarf satellites, 
we can include them in the BTFR. Fig. Q] shows the 
BTFR with the Local Gr oup dwarfs together with bright 
spirals (|McGaughl 120051) and gas rich, lat e type disks 
(|Stark et al.l 120091; ITrachternach et al.ll2009l) . The latter 
provide a calibratio n of the BTFR t hat is largely inde- 
pendent of the IMF (jStark et al.ll2009l ). The best fit slope 
is not meaningfully different from 4, so fixing the slope 
to this valu e and comb ining the gas rich galaxy d ata of 
iStark et al.l (|2009[ ) and ITrachternach et al.l (|20()9; leads 
to 

M b = AV C 4 (2) 



with A = 45 ± 10 M km _4 s 4 (jMcGaugh et al.l l20Tot) . 
Fig. Q] also shows the residuals after dividing out equa- 
tion $2$. The residual F b = M h j{AV 4 ) is the fraction of 
baryons a galaxy has relative to the expectation of the 
BTFR. 

The Local Group dwarfs diverge systematically from 
the BTFR. Figs. [H and [3] shows their residuals as func- 
tions of various quantities. Correlations are apparent 
with both luminosity and ellipticity, with fainter and less 
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Table 2 



Dwarf Tidal Radii and Ellipticities 


Dwarf 


ri/2 


r t,phot r t,D r t,M 


7 


c 


Fiefs. 



Carina 




,; t7 


846 ± 


106 


Draco 


291 ± 


14 


997 ± 


13 


Fornax 


944 ± 


53 


2854 ± 


161 


Leo I 


388 ± 


64 


850 ± 


63 


Leo II 


233 ± 


17 


556 ± 


28 


Sculptor 


375 ± 


54 


1758 ± 


115 


Sextans 


1019 ± 


62 


4000 ± 1250 


Ursa Minor 


588 ± 


58 


1495 ± 


171 


Bootes I 


321 ± 


28 


860 




Canes Ven. I 


736 ± 


17 


3000 




Canes Ven. II 


97 ± 


16 


500 




Coma Ber. 


100 ± 


13 


240 




Hercules 


304 ± 


26 


1500 




Leo IV 


151 ± 


39 


700 




Leo T 


152 ± 


21 


568 ± 


118 


Segue 1 


38 ± 


9 


160 




Ursa Major I 


415 ± 


60 


1400 




Ursa Major II 


183 ± 


33 


500 




Willman 1 


33 ± 


8 


110 




And I 


900 ± 


75 


2530 ± 


230 


And II 


1659 ± 


53 


4170 ± 


200 


And III 


638 ± 


77 


2120 ± 


360 


And VII 


1050 ± 


60 


4600 ± 


430 


And X 


448 ± 


8 


1480 ± 


150 


And XIV 


461 ± 


155 


1010 ± 


200 


References. 


— Deprojected half lig 


it an 



4220 ± 
4520 ± 
10350 ± 
10230 ± 
6170 ± 
4660 ± 
5900 ± 
5640 ± 
3920 ± 
10400 ± 
3060 ± 
1300 ± 
4240 ± 
2870 ± 
9470 ± 
630 ± 
4940 ± 
1730 ± 
780 ± 
25500 ± 
22300 ± 
13300 ± 
25700 ± 
10000 ± 
14500 ± 



380 
360 
720 
1680 
540 
510 
390 
490 
710 
770 
570 
220 
670 
1140 
1810 
150 
690 
400 
210 
2400 
1800 
3600 
3200 
2100 
4200 



1230 ± 
840 ± 
5890 ± 
6660 ± 
3400 ± 
1940 ± 
1270 ± 
890 ± 
330 ± 
2350 ± 
510 ± 
120 ± 
480 ± 
530 ± 
6180 ± 
38 ± 
420 ± 
103 ± 
79 ± 
17200 ± 
19200 ± 
10500 ± 
27900 ± 
5200 ± 
6900 ± 



250 
150 
1160 
1710 
610 
360 
230 
180 
60 
370 
150 
30 
120 
150 
1000 
12 
70 
28 
28 
1600 
1900 
1200 
3300 
800 
1300 



± 5 
l± 2.9 
± 10 
± 62 
±33 
± 8 



10 

6.9: 
22 
100 
78 
17 

2.0: 

2.6: 

1.5: 

8.0: 
17 

1.9: 

2.1 

9 

137 ± 23 
1.4 ± 0.7 
1.4 ± 0.7 
0.6 ± 0.4 
5.3 ± 4.6 



0.9 
1.3 
0.6 
3.0 
11 
1.1 
1.6 
7 



118 
56 
94 

193 
56 
81 



± 40 
±20 
±43 
± 83 
±32 
± 64 



0.33 ±0.05 
0.31 ±0.02 
0.31 ±0.03 
0.21 ± 0.03 
0.13 ±0.05 
0.32 ±0.03 
0.35 ±0.05 
0.56 ±0.05 
0.39 ±0.06 
0.39 ±0.03 
52 + 10 
0.36 ±0.04 
0.67 ± 0.03 

S 0.15 
0.29 ±S£ 

o.48 

0.80 ±0.04 
63 + 03 
0.47 ± 0.08 
0.22 ±0.04 
0.20 ±0.08 
0.52 ±0.02 
0.13 ±0.04 

0.31 ± 0.09 



1,6,11 
1,7,11 
1,6,11 
1,6,11 
1,6,11 
1,6,11 
1,6,11 
1,6,11 
1,7,11 
1,7,11 
1,7,11 
1,8,11 
1,9,11 
1,10,11 
1,7,11 
1,7,11 
1,7,11 
1,7,11 
1,7,11 
2,3,11 
2,3,11 
2,3,11 
2,3,11 
2,4,11 
2,5,11 



d p hotom etr ic tidal radii: 1 . IWolf et alj 1 120101 1 2. IKalirai et~ 



(2010 



19! IS 



3. IMcConnachie fc Irwrnl 02009) 4. IZucker et"aT1 1 1200711 5. IMaiewaki et alj 120071). E ll ipticity data: 6. IMatec 
7. iMartin et alj 1 1200831) 8. IMunoz et alj 120101) 9 jSand et alj l !2009al ) 10. ISand et all Il2009bf) ■ Tidal radii with dark 
matter and MOND and 7: 11. this work. 

Note. — All radii arc in parsecs. Columns: (1) Name of the dwarf satellite. (2) The deprojected 3D half light radius. 
(3) The photometric tidal radius. Cases lacking error bars arc approximate estimates only. (4) The tidal radius with 
dark matter from equation with rn — M^^. (5) The tidal radius in MOND from equation |5j with m — M^. (6) 
The number of orbits a star at the deprojected half light radius completes for every orbit of the dwarf about the host 
(equation[6j. (7) The ellipticity (e — 1 — b/a) of the dwarf as projected on the sky. (8) The source of the data. 



round dwarfs lying farther from the BTFR established 
for disk galaxies. The residuals also correlate with size, 
surface brightness, and metallicity. In these cases, the 
correlation is not as strong as with luminosity, and ap- 
pears simply to follow from the fact that these quantities 
correlate with luminosity themselves. The residuals cor- 
relate only weakly with current galactocentric distance, 
though there is a clear trend among the ultrafaint dwarfs 
of the Milky Way. This suggests a possible role for the 
o rbits of the dwar f s. 

IBellazzini et al.l ([1996D found a correlation of the sur- 
face brightnesses of dSph satellites with the combination 
of luminosity and absolute magnitude \My\ +6.41og(Z)). 
They argued that this suggested a role for environmental 
influences in shaping the structure of the dwarfs. We find 
that this combination of luminosity and absolute magni- 
tude also correlates with the degree of deviation form 
the BTFR (Fig. EJ). Indeed, it extends to include the 
ultrafaint dwarfs, which o f cours e were not know at the 
time that rBcllazzi ni et al.l (|1996t l identified this quantity 
a s relevant. 

IBellazzini et al.l (|1996l ) further define a dimensionless 
tidal force Ft,d as an indicator of the potential impor- 
tance of the gravity of the host galaxy on its satellites. 
For spherical masses in circular orbits, 



behavior. In equation ([3]), m is the mass of the satel- 
lite and r is its size while M is the mass of the host 
and D is the radius of the orbit. In Fig. [3] we use the 



deprojected 3D half light radius 



1/2 



as a measure of 



each dwarf's size, and set m = and Mh os t — V 2 D/G 
where V is the orbital velocity at the present location of 
the satellite. For the Milky Wa y we use the model ro- 
tation curve of McGaugh (2008) for the orbital velocity, 
which is typically ~ 190 kins" 1 at the d istance of most 
of the dwarfs (see also Kue et a l. 2008). For M31 we 
use Vm31 = 250 kms" 1 (|Chemin et alJl2009D . The larger 
Ft,d, the more susceptible a satellite is to the gravita- 
tional influence of its host. We find that the deviation 
from the BTFR correlates with Ft,d, a point we examine 
in greater detail in §§ 13.31 and !3.4l 



Ft,d — 



M fry 



(3) 



Other mass distributions follow the same dimensional 



3.1. Accuracy 

At first glance, the dwarfs follow the same trend as 
the spirals, just with larger scatter. Our first concern is 
therefore whether the data are accurate enough to fal- 
sify the null hypothesis that all of the dwarf Spheroidals 
are consistent with the BTFR. The classical dwarfs, for 
which the data are best, adhere fairly well to the BTFR. 
Indeed, Carina, Fornax, Leo I, Leo II, Sculptor, and Sex- 
tans follow the BTFR as well as rotating disks do. 

The M31 dwarfs also lie near to the BTFR, albeit with 
larger scatter. Three (And I, X, and XIV) are con- 
sistent with t he BT FR while the other three from the 
IKalirai et all (|2010l) sample (And II, III, and VII) sit 
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Figure 1. The masses and circular velocities ofgalaxies: the Baryonic Tully-Fisher Relation (left panel) and residuals from the BTFR 
(right panel) obtained by dividing out equation <J2J> (lines). Gray circl es are rotating disk galaxies. Dark gray po ints are star dominated 
spirals (McGaugh 2005) while light gray points are gas dominated disks (Stark ct al. 2009; Trachtcrnach ct al. 2009). The circular velocities 
of non-rotating galaxies are estimated as V c = V^^ios (equation Q}. Larg e blue squares are the classical dwarf Spheroidals and the small 
green squares are the ultrafaint dwarfs as tabulated by Wolf ct al. (2010). Leo T is the only dwarf spheroidal with gas; it is marked as a 
light blue star. Red triangles are the M31 dwarfs ijKalirai et aI.II20101 )~ 

somewhat above it. This might happen if equation (pj 
understates the circular velocity. The kinematic data we 
adopt for the M31 dwarfs are very recent. If we look 
back to a time when the data for the classical dwarfs 
were in a similarly early state, th eir scatter goes up 
(|Gerhard fc Spergell [19921: lMilgromlll995l) . We therefore 
consider the M31 satellites to be broadly consistent with 
the BTFR, given the uncertainties. And II is the most 
deviant case, being about 2a away. This is the situation 
for the most recent measurement, ai os = 7.3±0.8kms _1 
(jKalirai et al.ll2010[ ). This case illustrates the modest yet 
important way in which the data can evolve. If w e adop t 
(Tios = 9.3±2.7kms _1 as measured bv lCote et al.1 (|1999l ). 
no discrepancy with the BTFR is inferred. 

The majority of the MW ultrafaint dwarfs, together 
with Draco and Ursa Minor, lie systematically below the 
BTFR. Not all of the recently discovered dwarfs deviate. 
Canes Venatici I, Leo T, and Leo IV are all consistent 
with the BTFR within the errors. However, Bootes I, 
Canes Venatici II, Coma Berenices, Segue 1, Ursa Ma- 
jor I, Ursa Major II and Willman 1 all deviate substan- 
tially and significantly (> 2a) from the BTFR. This 
corresponds to deviation by at least an order of mag- 



nitude (F b ~ > 10). If a population of as-yet undiscov- 
ered ultrafaint dSphs exist with la rge dispersions and ex- 
tremely low surface brightnesses (|Bovill fc Ricottil 120091 : 
iBullock et al.l l2010). we expect them to deviate still far- 
ther from the BTFR. 

Hercules is an interesting case. It either deviates or it 
does not, depending on which publis hed velocity disper - 
sion we adopt. For consistency with iWolf et al.l ([2010), 
we use he re the velocity d i spersi on ai os = 5.1±0.9 kms -1 
found bv lSimon fc Gehal (I2007I) . For this aios, Hercules 
deviates substantially from the BTFR, being off by a 
factor iV 1 ~ 25 in mass, with a 68% confidence in- 
terval 7 < F b ~ < 45. Put this way, it would seem Segue 1, F b 
safe to say that it deviates by an order of magnitude. 
However, if instead we adopt the more recent measure- 



ment of oios = 3.7 ± 0.9 km s 1 reported by I Aden et al.l 
(|2009D . we find F" 1 « 7 with a 68% confidence inter- 
val 0.2 < Fj~ < 13. The central value again deviates 
by a large factor, and the uncertainties in the two mea- 
surements overlap. However, the ±1<t boundaries now 
encompass the BTFR, so we would no longer consider 
this to be a deviant case. 

The case of Hercules illustrates an important point. A 
small and formally consistent difference (1.4± 1.3 kms -1 
in this case) can have an impact on how we perceive 
the the results under consideration here. It is therefore 
important to take the formally significant deviations with 
a grain of salt. This is a rapidly evolving field; even 
with heroic efforts the velocity dispersion data for the 
ultrafaint dwarfs are not of the same quality as those 
for the classical dwarfs, nor can they be given the small 
number of stars in these systems. 

Given the cautionary tale of Hercules, one possible in- 
terpretation is that there are no real deviations from 
the BTFR. This hypothesis predicts that the scatter 
will decrease as the data improve, with the more de- 
viant cases migrating towards the BTFR. There is room 
to improve the d ata by weedi ng out nonmembe r stars 
(jSerra et al.ll2009f ) and binaries (jMinor et al.ll2010f ). both 
of which inflate ai os . Measuring proper motions to con- 
strain anis otropies will give a more accurate assessment 
of V c (e.g.. IStrigari et aLll2007|) . Of course, we must also 
consider the degree to which it is appropriate to compare 
ai os measured at small radii in dwarfs to V c measured at 
large radii in disk galaxies. 

There are substantial difficulties with attributing all 
deviations from the BTFR to errors. One is their sheer 
magnitude: the discrepancies Fb are an order of mag- 
nitude in some cases, and two orders of magnitude for 
some of the ultrafaint dwarfs. In the most extreme case, 



230. To place this object on the BTFR 



2 The difference in the velocity dispersions stems from how mem- 
bers of the system are identified, a systematic effec t that is no t 
represented by the formal uncertainty (see also Scrra ct al. 2009). 
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Figure 2. Residuals from the BTFR correlate with luminosity (top left), size (top right), surface brightness (bottom left), and metallicity 
(bottom right) in the sense that the faintest, smallest, dimmest, and most metal poor dwarfs deviate farthest from the BTFR. Symbols as 
per Fig. Q] 



would require that its baryonic mass increase by this fac- 
tor. Alternatively, since the BTFR is steep, a reduction 
in the observed a = 4.3 ± 1.1 kms" 1 (jGeha et all 120091) 
to aios = 1.1 kms -1 would also suffice. While a 3cr 
change may be conceivable in one object, it would have 
to happen in numerous cases. Moreover, the residuals 
are asymmetric: many more systems have Fb < 1 than 
Fb > 1. This asymmetry is larger than indicated by 
the formal errors, though we caution that the errors are 
asymmetric in the same sense. As always, a systematic 
error could be responsible, such as the inflation of o~i os 
by binaries. Any such systematic would have to similarly 
affect many different objects observed by a number of in- 
dependent groups while not affecting the objects that do 
adhere to the BTFR. 

Another problem with attributing deviations from the 
BTFR entirely to uncertainties is that the residuals cor- 
relate with the physical properties of the dwarfs (Figs. [2] 
and [3]). It seems unlikely that the observed correlations 
are mere flukes. If not, they will persist as data continue 
to accumulate. 

3.2. Gas Removal 



If the BTFR is a property shared by isolated galaxies 
when initially formed, we might seek to reconcile the Lo- 
cal Group satellites with it by extracting baryons before 
they form stars. Mechanisms to do this come in a variety 
of flavors. In this section we consider astrophysical effects 
that systematically suppress the cooling or retention of 
gas in dwarfs. 

3.2.1. Cosmic Reiomzation 

The reionization of the universe at z > 6 is one mech- 
anism by which the formation of stars in low mass halos 
might be suppressed. Reionization heats the gas, inhibit- 
ing further star formation in small h alos (V c < 20 kms" 1 
at the time of reinoization: see, e.g., IBullock et afll200a 
Kravtsov 2010). Halos that remain this small may never 
experience further star formation, pe rsisting as "fossil" 
dwarfs containing only anci ent stars (|Ricotti fc Gnedinl 
2005|; iBovill fc Ricottil 120091 ) . 

The results of sim ulations (jCrain et al.l 120071 : 
iHoeft fc Gottloeb"erll2010D are compared to the data in 
Fig. SJ Here we plot the detected baryon fraction, which 
is the fraction of baryons fd = Mb/ (fbM^oo) that are 
detected in any given halo relative to the cosmic baryon 
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Figure 3. Residuals from the BTFR correlate with indicators of environmental and and tidal influence such as distance from the host 
ga laxy (top left), shape ( as measured by the projected cllipticity of the dwarf, top right), the combination of luminosity and distance found 
by Bcllazzini ct al. (1996, bottom left), and a measure of the dwarfs' susceptibility to tidal influences (bottom right; equation [3} . The most 
deviant objects tend to be the closest (especially among the ultrafaint dwarfs) and least spherical. Leo T falls off the right edge of the top 
left plot with D = 407 kpc. Symbols as per Fig. [T] 

suggests a common origin. It appears that physics be- 
yond rcionization is needed to explain the observed trend 
in detected baryon fraction fd with halo mass. 

Reionization might act in addition to whatever physics 
drives the main trend in Fig. |4] by causing a sharp low- 
velocitj{3 cut-off. This might explain the residual corre- 
lation with luminosity exhibited by the ultrafaint dwarfs 
(at least those with Ly < 10 5 L Q in Fig. [5]). However, it 
does nothing to explain the trend of the residuals with 
ellipticity (Fig. [3]). 



fraction f b (|McGaugh et al.ll2010D . Note that f d differs 
from the BTFR residual F b . 

R eionization models like tho se of iCrain et al.l ()2007f ) 
and iHoeft fc Gottloeberl (|2010f l result in a sharp trun- 
cation of the cold baryon content of halos at a particu- 
lar mass scale. Above this mass scale, reionization does 
nothing to prevent the baryons from cooling and forming 
stars. In contrast, the data deviate smoothly and grad- 
ually from the universal baryon fraction, not abruptly 
as suggested by models consisting of reioni zation only. 
Thus, while the model of ICrain et al.l (|2007t) appears to 
work well for V c < 20 kms" 1 , it does not explain dwarf 
Irregulars with V c > 30 kms" 1 . At this slightly larger 
scale, reionization is not expected to impede the conden- 
sation of cool gas and subsequent star formation, yet the 
mass of baryons detected in the form of stars and cold 
gas is an order of magnitude shy of the cosmic baryon 
fraction. 

We might suppose that there is some fundamental dif- 
ference between dSph and dlrr galaxies that separately 
determines their location in Fig.|4] However, the continu- 
ity of the BTFR to most of the classical dwarf spheroidals 



3.2.2. Stellar Feedback 

Stellar feedback is frequently invoked as a mechanism 
to remove gas from galaxies or inhibit its accretion in 
the first place. Intense radiation fro m young stars (per- 
haps Pop. Ill stars at early times, e.g.. lRicotti fc Ostrikerl 
2004) may heat the surrounding gas and prevent its infall. 
Winds driven by supernovae following intense bursts of 

3 An interesting test of the idea that reionization imposes a sharp 
cut-off to the existence of low mass galaxies would be if there exist 
any galaxies of exceedingly low rotation velocity (V c < 5 kms -1 ) 
that persist in following the BTFR, especially if they were gas rich. 
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Figur e 4. The dete cted baryon fractions = Mh/if^Msoo) 
(McGaugh ct al. 2010) of galaxies (symbols as per Fig. [lj fall well 
below the cosmic baryon fraction /(, (dotted line). The ultrafaint 
dwarfs appear to follow a steeper relatio n than the BTFR (thin 
solid line): the dashed line is the fit from McGaugh et al. 120101) 
before the M31 data (Kalirai ct al. 2010) were available. The 
heavy solid lines illus trate the expected effects of cosmic reioniza- 
tion (descending line: Crain ct al. 2007; Hocft & Gottlocbcr 2010, 
S-shaped line). This alone cannot simultaneously explain both 
dwarfs and disks. 



star formation may impart sufficient mechanical energy 
into the interstellar medium to sweep up and expel sub- 
stantial amounts of g as (e.g.. lMac Low fc Ferraralfl999at 
IGnedin &T Zhao 2002)). In either case, the prodigious en- 
ergy produced by massive, short-lived stars is invoked to 

alter the global properties of ga laxies. 

These ideas go back to at least iDekel fc SilkJ (If 9861) and 
have been invoked many times (e.g.JMac Low fc Ferraral 
f 999bt iNataraianllf 999[ lEfstathioull2000t Ivan den Boschl 
2000t iMaver fc Moorel 120041: iRicotti fc Gnedinl 1 2005 ; 
Mashchenko et all 120081; iDutton & van den Boschl 12009 ; 



Governato et al.ll20f Ct iSawala et al.l 120101 ) . One appeal 



ing aspect of this picture is that the energy injected by 
star formation events depends on the local star formation 
efficiency, while the ability of a halo to retain baryons de- 
pends on the depth of its potential well. It is therefore 
natural to expect a trend like that seen in Fig.2]with the 
retained baryon fraction increasing with halo mass. Star 
formation happens in a local way, but global baryon re- 
tention becomes increasingly difficult in smaller galaxies. 
When the potential well becomes very shallow, as is the 
case for dwarf spheroidals, stochastic effects may come 
to dominate completely such that Fb varies widely at a 
given V c . 

While qualitatively attractive, quantitative implemen- 
tation of feed back in realistic simul ations is notoriously 
difficult (e.g.. lGovernato et al.ll2010h . The gas physics is 
uncertain, and even the depth of the potential well is un- 
clear. In the ACDM structure formation paradigm, large 
galaxies are constructed from the mergers of smaller ones. 
The trend of fd with potential well depth only follows 
naturally if most of the star formation activity responsi- 
ble for the feedback takes place after the establishment 
of the bulk of the potential well depth as indicated by the 
circular velocity measured now. If instead much of the 
action occurs in small protogalactic clumps before the 
establishment of the large halos we imagine to retain the 



gas, it becomes less obvious that the ultimate potential 
well depth is relevant. Indeed, as our nearest examples 
of the fossils of protogalactic clumps, the large scatter 
in luminosity and fd at a given mass exemplified by the 
dwarf spheroidals more nearly represents the stochastic 
mess that we would naturally expect from this process. 

Given that stellar feedback is a stochastic process, any 
relation that it establishes will likely have a great deal of 
scatter. In the context of the Local Group dwarfs, this 
is quite reasonable: there is a lot of scatter about the 
BTFR (Fig. [T]). This does not persist to larger masses 
(V c > 30 kms" 1 ). For the gas dominated late type disks 
popu lating Fig. H (jStark et all 120091 : iTrachternach et al.l 
120091 ) . the data are consistent with a relation with nearly 
zero intrinsic width. These data already place an uncom- 
fortably tight limit on the t riaxiality of dark matter halo s 
(|Trachternach et al.l 12001 iKuzio de Narav et al.l [2009h : 
there is not much room for additional stochastic effects. 

Note also that rotating galaxie s adhere to the BTFR 
regar dless of their gas fraction (jAnderson fc Bregmanl 
120101 ). This is contrary to the natural expectation of 
feedback scenarios in which galaxies that have converted 
more gas into stars have experienced more star formation 
and hence more feedback. For the dwarf spheroidals, the 
increased scatter about the BTFR could be the signature 
of a stochastic effect, if not for the fact that the residu- 
als are not random as they correl ate with other physi cal 
properties (Figs. [2]and|31 see also Walke r et al.l 1 201 01). 

Recent simulati ons (e . g .. iStinson et al.l [20091 : 

iSalvadori fc Ferraral 1200a ISawala et al.l 120101) of 
low mass galaxies appear to do a reasonable job of 
reproducing some of their observed properties. It will 
be interesting to see if they are equally successful in 
matching the correlations of residuals from the BTFR. 
The correlations with luminosity and metallicity may 
indeed follow, as the chaos of feedback can easily lead 
to a wide range of luminosity at a given halo mass, and 
supernova winds may carry away much of the metals 
they produce. The correlation with ellipticity and tidal 
susceptibility is less obviously explicable in this way. 

3.2.3. Ram Pressure Stripping 

Another mechanism capable of removing cold gas be- 
fore it can form stars is ram pressure stripping. In this 
hypothesis, the initial condition is gas rich galaxies that 
adhere to the BTFR. As these objects fall into larger 
structures like the halo of the Milky Way, hot gas in the 
host's halo ablates the cold gas of the infalling satellites. 
The resulting object is under-luminous with respect to 
the BTFR because many baryons were removed before 
they could form stars. 

The Magellanic str eam may provide an examp le of this 
process m progress (jMastropietro et al.l 120091). though 
tidal effects might also suffice (jConnors et al.1 120061 : 
iStanimirovic et ail 120081) . Certainly it is striking that 
all of the nearby ( D < 250 kpc) dwarf sphe roidals are 
devoid of gas (e.g.. lGrcevich fc Putm an 2009). Provided 
that there is sufficient hot gas in the halo, ram pressure 
stripping should occur. 

In this hypothesis, we associate the stripping of the 
cold gas baryons with the epoch when the object fell into 
the larger structure. The galaxies that fell in first would 
have had the least opportunity to form stars and would 
have lost the most gas. This has two potentially testable 



consequences. First, the amount by which a dwarf devi- 
ates from the BTFR is an indicator of the infall epoch. 
The earlier it fell in, the sooner gas was stripped and 
star formation truncated. Secondly, because of the trun- 
cation of star formation, the age of the stars should also 
correlate with the factor by which the dwarf deviates. 
Dwarfs that fell in early should contain only old stars. A 
natural corollary is that the stars in the first dwarfs to 
fall in should be low metallicity. Indeed, since these ob- 
jects would have had little time to self-enrich, we would 
expect the [a/Fe] ratios to be enhanced while [Fe/H] is 
low. We would thus expect [a/Fe] to correlate with the 
residuals F^. 

The degree of deviance F}, should reflect the order of 
infall. The dwarfs that deviate most presumably do so 
because they fell in earliest. If this holds, then Segue 1 
and Ursa Major II fell in first, followed by Willman 1, 
Bootes I, Ursa Major I, and Coma Berenices. Hercules, 
Canes Venatici, Draco, and Ursa Minor are relative new- 
comers, while Leo T, which is the only MW dSph with 
detectable HI gas, has yet to be stripped. 

The prediction that deviance correlates with infall time 
is not perfect, as it depends on the star formation history 
of the dwarf prior to infall. An object that converted 
most of its cold gas into stars quickly before infall would 
not deviate much from the BTFR even if the infall were 
relatively early. However, no star formation can persist 
after ram pressure stripping is complete. We do therefore 
expect the dwarfs composed of the oldest stars to be the 
first to have fallen in. 

An interesting corollary has to do with the epoch of 
most recent star formation. Some of the c lassical dwarf s 
have had complex star forming histories (jGrebell [l997f l. 
and are not composed solely of ancient stars. In some 
cases , the last epis ode of star formation is rather recent 
/e.g. lGrebellll997t Hernandez et al1l2000l: lDolphirJl2002t 
iHelmil I2008t iMartin et all I2008U iTolstov et al.l 12009) . 
This has been a puzzle, since, with the exception of Leo 
T, none of these objects contain substantial quantities of 
cold gas at present. The truncation of star formation by 
infall predicted by the ram pressure stripping hypothe- 
sis provides a natural explanation for this puzzle. The 
dwarfs with recent star formation episodes are merely the 
most recent to have fallen in and had their gas stripped 
away. 

For this process to be generic, the halos of host 
galaxies need to c o ntain a good deal of hot gas. 
iGrcevich fc Putmanl (|2009l ) infer a density of hot gas 
n > 2 x 10~ 4 cm~ 3 out to at least 70 kpc around 
the Milky Way in order for ram pressure stripping to 
have effectively ablated t he local population of dwar f 
spheroidals. In contrast, lAnderson fc Bregmanl ([2010D 
place a limit n < 7 x 10~ 5 cm~ 3 at 50 kpc based on 
pulsar dispersion measures towards the LMC and other 
methods. This would not be enough for ram pressure 
stripping to have played a major role in shaping the his- 
tory of satellite dwarfs. The mass and extent of the hot 
gas halos around spiral galaxies is an important but ill- 
constrained piece of the puzzle. 

The effects of ram pressure stripping are not limited to 
the small scales of interest here. It is obse rved to affect 
spira l galaxies falling into rich clusters (e.g.. lChung et al.l 
2009). If this mechanism causes dwarf spheroidals to 
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Figure 5. The ellipticity of Local Group dwarf Spheroidals as a 
function of distance from their host galaxy. Symb ols as per Fig. [l] 
The open square is Sagittarius (Ibata et al. 1997). Leo T (e = 0.29 
at D = 407 kpc) falls off the right edge of the plot. 

deviate from the BTFR in the Local Group, we might 
expect it to have the same e ffect on spirals in clusters . 
However, this is no t observed ()Guhathakurta et al.lll98cl 
IChung et afll2009h . Since hot gas is definitely present in 
rich clusters, and cold gas is observed to be stripped from 
spirals there, one might expect the effect to be more pro- 
nounced in clusters rather than less. However, stellar 
mass dominates the baryon content of the infalling spi- 
rals, which typically have modest gas fractions (~ 0.1 - 
0.2). This might suffice to mask the effect, though an 
investigation of the impact of different environments on 
the BTFR would clearly be interesting. 

Ram pressure stripping potentially provides a good 
way to strip gas at an appropriate time to allow for recent 
star formation in some dwarf spheroidals that possess no 
cold gas now. It would naturally account for the correla- 
tion of BTFR residuals with luminosity and metallicity 
if enough gas is stripped. It does not, by itself, provide 
an obvious explanation for the correlation with ellipticity 
and tidal susceptibility. These are features of the stellar 
distribution that are difficult to explain with any mech- 
anism that acts solely on the gas. 

3.3. Stellar Stripping 

Tides are a mechanism that acts on the stars directly. 
The correlation of BTFR residuals with ellipticity and 
tidal susceptibility as well as with luminosity would oc- 
cur naturally if some of a dwarf's original stars have 
been stripped away by tidal perturbations as it orbits the 
Milky Way. Under this hypothesis, the deviant dwarfs 
originally adhered to the BTFR, but have had much of 
their original luminosity removed. Much of the stel- 
lar halo of the Milky Way might hav e been built up 
this way through hierarch ical merging () Johnston! 119981 : 
IBullock fc Johnstonl[2005l) . 

Analyses of the dynamics of the dwarf spheroidals typ- 
ically assume spherical symmetry. The large ellipticitics 
of many of the dwarfs call this assumption into ques- 
tion. It is manifestly untrue in Canes Venatici II, Ursa 
Minor, Ursa Major I and II, and Hercules, all of which 
have e > 0.5. In this context, it is important to deter- 
mine whether some dwarfs show dynamical evidence of 
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disruption fe.g.. IMunoz et al.ll2008l 12010ft . 

The fact that the degree of deviation from the BTFR 
correlates with ellipticity (Fig. [3]) may be a clue that tidal 
disruption is playing a role. A further clue is provided by 
Fig. [3 which shows that the ellipticity of the dwarfs also 
correlates with distance from their host. While we do 
not expect all dwarf Spheroidals to be perfectly round, 
neither do we expect their intrinsic shape to depend on 
distance unle ss some environment al factor is playing a 
role (see also iBcllazzin i et al.lll996h . 

If ellipticity is an indicator of tidal disruption then the 
correlati on of deviations with luminosity follows natu- 
rally (cf. IPiatek fc Prvorl[l995t IMunoz et al.ll2008ft . The 
more disrupted a galaxy, the more elliptical it becomes, 
and the more stars it has lost. The more stars a 
dwarf has lost, the further it deviates from the BTFR. 
An obvious observational test of this hypothesis would 
be to search for streams associ a ted with the deviant 
dwarfs (e.g.. IMunoz et al.l 120051: iFellhauer et all 120071: 
iSales et all 120081: iNewberg et al.l 12010ft . The age and 
metallicity of stars in these streams should be consistent 
with having been drawn from the distribution present in 
the parent body. 

Sagittarius is one example where this has clearly hap- 
pened. It is now out of equilibrium, so we would not 
expect it to reside on t he BTFR. If we ad opt a lumi- 
nosity of - 2 x 10 7 L (|Mateo et al.l 119960 . the BTFR 
predicts ai os ~ 16 km s" 1 , somewhat larger than the ob- 
serve d ~ 11 kms -1 (|Ibata et al.l Il997t iBellazzini et al.l 
2008). Given the uncertainties, it is not obvious that 
this constitutes a significant discrepancy. However, if 
the larger total lumin o sity e stimate (~ 10 s L Q ) of 
iNiederste-Ostholt et al.1 ()2010h is correct, then Sagittar- 
ius should initially have had ai os ~ 24 kms -1 . We need 
to understand how luminosity and velocity dispersion 
evolve as disruption takes place. 

A satellite in orbit around a much lar ger parent galax y 
is obviously subject to tidal disruption | Johnston! 11998), 
but it is less obv ious how ai ns will evolve as mass is 
stripp ed (see, e.g.. lKlimentowski et al.l[2009l : lLokas et al.l 
[20nj. Naively, we would expect the initial stripping to 
primarily affect the outer dark halo first. Stars would 
not be stripped in significan t numbers until much o f the 
dark halo was already gone (|Peftarrubia et af1l2008[ ). 

We can be more quantitative about the propensity of a 
dwarf to be stripped by returning to the subject of tidal 
susceptibility (e quation EQ) . It is conventional to define 
the tidal radius ((KeenarfosTallbh as 
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This is approximately the radius where a star contained 
in a satellite of mass m is subject to being lost to the 
host mass M, though the details depend on the orbit. 
Equation Q implicitly assumes circular orbits; eccentric 
orbits are most susceptible to tidal stripping near peri- 
center where D — > a(l— e) (where a is the semi-major axis 
and e is the eccentricity of the orbit). Equation (0| pro- 
vides an estimate for when stripping may occur, bearing 
in mind that the current galactocentric distances of the 
dwarfs could be larger than their pericenter distances. 

Fig. [6] shows the sizes of the dwarfs in terms of tidal 
radii. For a star orbiting at the deprojected 3D half 



light radius, we set m = M x / 2 = 3of os ri/ 2 /G (Wo lf et all 
12010ft . The distance D is the galactocentric distance of 
each dwarf from its host, and the mass of the host is the 
mass contained within D: M = V 2 D/G For the Milky 
Way we use VkwP) from iMcGaughl (j2008) while for 
M31 we assume that a constant velocity suffices at the 
location of the M31 d warfs and adopt Vm31 = 250 kms -1 
(IChemin et alJl2009h . 

The tidal radii of the dwarfs are typically an order 
of magnitude or more larger than their half light radii. 
By this criterion, the stars of all the dwarfs are safely 
cocooned deep inside their dar k matter halos and ar e 
in no danger of being stripped (Pcnarru bia et al.l l2009). 
The vast majority of the dark matter must be tidally 
stripped before the sta rs become exposed to stripping 
(|Penarrubia et aLl l2008). It therefore appears that these 
dwarfs should not currently be losing stars. This is some- 
what surprising, given the large deficit Fj, of some of these 
dwarfs, and the correlation with tidal susceptibility seen 
in Fig. H 

Photometric tidal radii are less certain than the depro- 
jected 3D half-light radius plotted in Fig. [6] but correlate 
very well with it. On average, r t . p hot = (3.1 ± 0. 8)7-1/2 
for the dwarfs with photometrically well measured tidal 
radii. We mark this location in Fig. [5] The mean size 
of the dwarfs is only 7*1/2 = 0.06r tj £>. To bring the pho- 
tometric tidal radii into agreement with the dynamical 
tidal radii computed with equation Q would require an 
adjustment to the mass of nearly two orders of magni- 
tude: m « O.OIM1/2. This seems like a lot to ask, but 
there are several possibilities. 

First, M1/2 ma y be overestimated because the dis- 
rupting systems are not in equilibrium. The true mass 
might be less, allowing tidal escape. In this case the 
velocity dispersion is artificially inflated by escaping 
stars. Though qualitatively attractive, this effect is prob- 
ably not la r ge en ough. From numerical simulations, 
lLokas et al.l ( 20101) estimate that the inflation of the ve- 
locity dispersion due to stripping should result in an over- 
estimate in mass of at most 60%. Nevertheless, one can 
imagine that the amount of dark matter is less than it 
might seem. 

In this context, it is worth noting that at least some 
dark matter is required. Equation (|4]) differs from equa- 
tion (J3j) only by a factor of 3 1 ' 3 , yet the result in the left 
panel of Fig. [6] differs from that in the lower right panel 
of Fig. [3] because there we set m = M5. Dwarfs lacking 
in dark matte r are quite susceptible to tidal influences 
(|Kroupalll99"7l) . Indeed, using the baryonic mass Mb in 
place of the dynamical mass M x / 2 in equation ((4]) results 
in tidal radii that are much smaller than the observed 
half light radii. If this were the case, the dwarfs should 
have disintegrated long ago. 

A second possibility is that Mi/ 2 is correct and the 
observed dwarfs are not yet the parents of any streams. 
The streams observed in the halo would have been pro- 
duced by the dissolution of other dwarfs. These objects 
dissolved in the past and are now completely gone. While 
this may have happened, it docs nothing to explain the 
deviations Ff,, nor their correlation with ellipticity. 

A third possibility is that the orbits of the dwarfs 
are highly radial, a situ ation that is wide l y ex- 
pected in ACDM (e.g., iBullock fc Johnston! [20051 : 



10 




o 
o 



CO 

I 

o 



I 

o 




10" 



10" 



10 u 



1/2 



At. 



M 



Figure 6. The residuals from the BTFR (symbols as per Fig. [TJ plotted against the ratio of half-mass radius to tidal radius in the case 
of dark matter (left) and MOND (right). The plots are scaled identically. The tidal radius depends on the masses of satellite and host 
(equations [4] and [5} . Masses are computed dynamically for the case of dark matter (left: m = M1/2) and using only baryonic mass in the 
case of MOND (right: m = Mi). The arrow marks the location where the average photometric tidal radius equals the computed tidal 
radius. In the case of MOND, the location where dwarfs deviate form the BTFR corresponds well to the observed photometric tidal radius, 
and the amount of deviation correlates with the size of a dwarf relative to its MONDian tidal radius. 



Klimentowski et al l 120091: IPenarrubia et al.l 120091 : 
Lokas et all 120101: ISales et all I2010D . If we replace D 
with a(l — e), rt.D shrinks as e — » 1. To explain Fig. 
we need the deviant dwarfs to be on orbits with e > 0.9. 

Most of the damage would occur during pericenter pas- 
sage, when the least bound stars are lost. The typical 
dwarf is not currently close to its pericenter, but may 
bear the scars of its close passages in the form of reduced 
luminosity and deviation from the BTFR. In this case we 
imagine that the streams in the halo have been built up 
by stars lost during successive pericenter passages. 

This last possibility provides an appealing explanation 
for the correlation of F\, with ellipticity as well as lu- 
minosity. However, we caution that the velocity disper- 
sion is subject to evolution during this process as well as 
th e luminosity and shape o f a dwarf. In the simulations 
of IPenarrubia et"al . (2008), the velocity dispersion de- 



creases as luminosity is lost, in a proportion that almost 
maintains the slope of the BTFR. Such an evolution will 
not reproduce the trend observed here, which requires 
a more rapid loss of luminosity with respect to veloc- 
ity dispersion. This motivates further simulation work 
that probes a wider range of initial conditions and halo 
models. 

In the context of our empirically motivated stellar 
stripping hypothesis, we expect three broad classes of 
objects. Dwarf satellites that have not yet ventured too 
close to their massive host will reside on the BTFR. The 
classical dwarfs (excluding Ursa Minor and Draco) and 
the dwarfs of M310 fall in this category. Dwarfs that de- 
viate from the BTFR do so because they have lost some 
stars during pericenter passages. Greater deviation is a 
sign of greater damage, from either closer pericenters or 
multiple pericenter passages. Ursa Minor, Draco, and 
most of the MW ultrafaint dwarfs fall in this category. 
Finally, some of the observed streams may have come 
from parent bodies that have been totally disrupted by 
this process. 

4 Examination of the properties of the dark matter halos of the 
M31 dwarfs suggests the opposi te, that they may have been more 
subject to tidal evolution (Walker ct al. 201(31). 



The stellar stripping hypothesis makes several predic- 
tions. The dwarfs that deviate from the BTFR should 
be the parents of stellar streams. These stars should 
have properties consistent with the parent body: their 
ages and abundance patters should be consistent with 
the same distributions. The streams should be in orbits 
identifiable with their parents. It appears likely that 
stars would largely be liberated during pericenter pas- 
sage, so there is in principle a further prediction about 
the phase at which stars were injected into the stream. 
Pericenter passage is a brief portion of the orbit, so we 
would not expect any of the observed dwarfs to currently 
be in the process of losing stars. At their present galacto- 
centric distances, all of the dwarfs should be safely bound 
in what remains of their dark matter cocoons. Observa- 
tion of streams being liberated at present would require 
some further interpretation (§ 13 .4|) . 

In order to liberate stars as envisioned under the stellar 
stripping hypothesis, the disrupting satellites need to be 
on highly elliptical orbits. Though we do not yet know 
the orbits of all the dwarfs, there are s ome constraints. 
Carina has e ~ 0.67 dPiatek et ail 20031 ) and Ursa Minor 
has e w 0.39 (|Piatek et al.l 120051 ) . Carina adheres to 
the BTFR but Ursa Minor does not, so their observed 
orbits are the opposite of what we would expect in this 
hypothesis. While disfavored, a pericenter passage that 
is sufficiently small (~ 20 kpc) to liberate some stars 
from Ursa Minor is not yet excluded by the orbital data 
(|Lux et alJ[2010h . Even if this were the case it would 
remain curious that Ursa Minor deviates from the BTFR 
and Carina does not. Clearly, better orbital constraints 
are desirable, as are searches for liberated stars. 

3.4. MOND 

An alternative to the ACDM paradigm is the Modified 
Newto nian Dynamics (MOND) hypothesized bv lMilgroml 
(1983). In MOND, there is no dark matter; rather, 
the apparent need for dark matter stems from a change 
to the force law that occurs at an acceleration scale 
ao < 1 A s~ 2 . Above this scale, all is normal: the 
gravitational acceleration can be calculated from the ob- 
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served distribution of baryonic matter and Newton's Law 
of Gravity: g = 9n- Below this scale (a <C ao), the 
modification g — > ^Ja$gN applies. This idea has its prob- 
lems ([Clowe et al.ll2004lAngus et al.ll2008D . but has also 
had more success than seems to be widely a ppreciated 
(jSanders fc McGaughlEfflOl : IBekensteinll2006l ). 

Dwarf spheroidals provide a strong test of MOND be- 
cause their low surface densities place them deep in the 
modified regime. In this regime, the BTFR is an abso- 
lute consequence of the force law for isolated objects: 
g = y/a^g~N leads to a^GMb = for point masses. 
The gas dominate d dwa rf Irregular galaxies studie d by 
iStark et al.l ((2009) and iTrachternach et al.1 (|2009f ) fall 
along the BTFR, consistent with MOND. Indeed, this 
is confirmation of a prediction that has zero free param- 
eters: the gas mass dominates Mb and these galaxies fall 
along the MOND prediction, which is indistinguishable 
from the best fit BTFR. 

The deviation of the Local Group dwarf satel- 
lites fro m the BTFR would a p pear to falsify 
MOND dGerhard fc Spergel [l99l IMilgroml 
Sanchcz-Salc edo fc Hernan dezl 120071: lAngusI 



1995 



2008 



Serraet al.1 120091: IHernandez et all 120101) . However 
the Local Group dwarf satellites are not isolated. In 
MOND, the criterion for isolation is whether the internal 
acceleration of an object, its self- gravity gi n t, exceeds 
that imposed by external systems, g ex t- If gint < 9ext, 
the external field modifies what would otherwise happen 
to the same object in isolation. This external field effect 
is a unique feature of MOND that has no analog in the 
conventional context. 

Among the dwarfs considered here, Leo T is the only 
example that is isolated enough to be in the pure MOND 
regime such that g ex t < gint < ^o- As such, it 
should obey the BTFR, and within the errors, it does. 
Sagittarius is an example that is clearly not isolated, 
with gi n t < g ex t < Qn- This pl aces it in the quasi- 
Newtonian regime ( Milgrom 11986!) where the mass es- 
timator is different (Milgroml ll995! ) and the BTFR is no 
longer absolute. The quasi-Newtonian mass estimator 
basically just corrects the normal Newtonian mass es- 
timator by a factor G — » G e // = Gao/g ext (|Milgroml 
119861) . Using a = 1.2 A s" 2 ([Begeman et al.lll99l!) and 
_ t/2 /n „,uv, y MW = 210 kms" 1 (jXue et al.l 
at galactocentric D = 19 kpc 
the mass-to-light ratio of Sagit- 
1.2 M fT> /L^ for the total luminosity esti- 
mated by |!EatairaDIl997). This is surprisingly reason- 
able, but would be reduced if we adopted th e large r lumi- 
nosity estimate of Niedcrste-Ostholt et al. (2010). Per- 
haps this reasonable seeming mass-to-light ratio is a fluke 
since all the same assumptions apply to the analysis here 
as in the Newtonian case: stability and sphericity. Nei- 
ther condition is satisfied in Sagittarius, which is highly 
elliptical and far from equilibrium. 

For the classical dwarfs excluding Draco and Ursa 
Minor, the inferred ma ss-to-light ra tios are reasonable 
for stellar populations (LAngugJ [200l ISerra et al.l 120091) . 
IHernandez et al.l (|2010f ) even find that lower mass-to- 
light ratios tend to occur in dwarfs with more recent 
star forming events, as expected from a stellar popula- 
tions perspective. This is not trivial, as the mass-to- 
light ratio computed in MOND is extremely sensitive to 
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the accuracy of the data simply because the mass de- 
pends on a high power of the velocity dispersion. Even 
small errors in identifying m embers can have a notice- 
able effect (|Serra et al.l 120091 ) . So far, however, Draco 
and Ursa Minor persist in being problem cases with 
uncomfortably high (~ 10 M^ /Lm) mas s-to-light ratios 
(IGerhard fc SnergeHfl^il IMilgroml 1 19911 ). 

For the ultrafaint dw arfs, the mass-to-light ratios 
are n ot at all reasonable (jSanchez-Salcedo fc Hernandez! 
I2007D . Though the BTFR is no longer absolute when the 
external field do minates, it shou ld remain a reasonable 
proxy for MOND (|Miigromll986!) . In this approximation, 
the inverse of the deviations F b ~ are proxy estimates for 
the mass- to-light ratios in MONE0 as they measure how 
re mote each dwarf is from the mean TY =_L3 M q /Lq 
of lMateo et"aH (|1998| ) and lMartin et all (|2008aft . From 
this perspective, Y]f « F^ 1 > 10 is quite unreasonable. 
The cases with Fb ~ 100 are simply absurd. 

The ultrafaint dwarfs may therefore be fatal to MOND, 
provided that both their data and the assumptions un- 
derlying the analysis are valid f§ 13. ip . The assumption 
of sphericity is obviously not satisfied in many of the 
most deviant cases, though it seems unlikely that geo- 
metric corrections alone could explain the large devia- 
tions seen in these objects. The assumption that the 
dwarfs are in dynamical equilibrium is perhaps more 
consequential. If they are, they should not exhibit dis- 
crepancies of one or two orders of magnitude. If they 
are not, then MOND may become similar to the stellar 
stripping hypothesis I3.3|) . though not identical. Tidal 
forc es are stronger in MOND than in conventional grav- 
ity ([Brada fc Milgromll2000b[ ) as the long range boost to 
the effective fo rce becomes more import ant than the ex- 
tra dark mass (| Angus fc McGaughll2008l) . This opens the 
possibility that stripping could be occurring at present. 

We can estimate the tidal radii of the dw arfs in MOND 
j ust a s we did conventionally in M3.31 IZhao fc TianI 
(|2006l) find that the expression for the tidal radius in 
MOND is nearly identical to that in the conventional 
case (equation!?]), differing only by a factor of (2/3) 1 / 3 : 



D (—) 

\2MJ 



1/3 



(5) 



What is really different here are the masses m and M 
of the satellite and host. Here we simply set m = Mb, 
the baryonic mass of each dwarf estimated from its lu- 
minosity as before. For the host mass M, we take 
M M w = 6 x 10 1 Mq for the Milky Way (|Flvnn et al.1 
[20061: IMcGaughl [2Q08T) . and M M31 = 2M MW for M31 
( Chemin et al.ll2009D . Note that these are the baryonic 
masses of the two host galaxies, as there is no dynamical 
dark matter in MOND. 

Contrary to the case for dark matter, the residuals 
from the BTFR correlate strongly with the size of the 
dwarfs as measured by MONDian tidal radii (Fig. [6]). 
The closer the half light radius is to the tidal radius, 
the further the dwarf deviates from the BTFR. Intrigu- 

5 The problem in these cases is lesse ned but not eliminat ed if we 
adopt the lower velocity dispersions of Walker et al. (2009). 

In this regard, we confirm the results of 



Sanchcz-Salcedo & Hernandez ( 2007) for the MOND mass- to-light 
ratios of these dwarfs. 
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ingly, the size scale at which the deviation occurs corre- 
sponds well to the photometrically measured tidal radius: 
ft,M ~ ft,phot- Galaxies that adhere to the BTFR are 
safely within their tidal radii such that r t , p h t < rt,M 
while the deviant cases typically have rt, p hot ^ ft,M- 
Dwarfs with a larger fraction of their stars exceeding the 
tidal radius appear to have lost more light. This implies 
that tidal disruption may be the cause of the deviance of 
Draco, Ursa Minor, and the MW ultrafaint dwarfs. The 
MONDian interpretation of these systems would thus ap- 
pear to be that they are in the process of being tidally dis- 
torted and ultimately shredded in the field of the Milky 
Way. The small F b are not large MONDian M*/Ly be- 
cause the assumption of stable equilibrium does not hold. 
Rather, a small Ft> is an indicator of the degree of tidal 
disruption. 

The deviations Fb in Fig. [6] appear to follow a straight 
line that continues from well below the BTFR to a bit 
above it. Satellites should not have Fj, > 1 as a result 
of tides in MOND. Rather, once objects become isolated 
enough, they should adhere to the BTFR with no net 
residuals. So this plot should go up to unity, then re- 
main there. Presumably this is what we would see if 
we could include the isolated disk galaxietQ from Fig. [TJ 
The perception of a linear rise to Ft, > 1 is due to three 
satellites of M31: And II, III, and VII. For these three 
objects, the implied MOND mass-to-light ratios are un- 
comfortably low: T» « 0.2 M Q /L Q would be required to 
bring them onto the BTFR. Given the uncertainties in 
both the analysis and the data (§ I3.1|) . we hesitate to put 
too much weight on these few cases. 

The correlation seen in the MONDian portion of Fig. |5] 
is quite strong. With the exception of the most extreme 
outlier, it is consistent with little intrinsic scatter. How 
great the intrinsic scatter is depends sensitively on what 
we take for the uncertainty in T*. We would not expect a 
perfect correlation, as there should be some variation due 
to the orbits of the dwarfs. The expression for the tidal 
radius implicitly assumes circular orbits; more generally 
we should replace D — > a(l — e). Dwarfs on radial orbits 
should be more susceptible to disruption, while retro- 
grade orbits may be somewhat more stable. While it is 
probably an over-interpretation of the data, we should 
at least point out that this effect might be perceptible. 
Two or three of the ultrafaint dwarfs (Canes Venatici II, 
Willman 1, and perhaps also the most discrepant case, 
Segue 1) sit below the main correlation. One might in- 
fer from this that they are on more eccentric orbits than 
the other dwarfs. Sextans, on the other hand, remains 
consistent with the BTFR in spite of having a half light 
radius approaching its tidal radius. This situation might 
persist longer if it happens to be on a retrograde orbit. 
Another possibility is that Sextans is currently making 
its first close approach, so has not y et suffered disruption. 
Indeed, iBrada fc MilgromI (|2000al ) show that dwarfs on 
mildly elliptical orbits may not disrupt at all if they are 
only exposed to pronounced tidal effects for a brief period 
near the pericenters of their orbits. 

The numerical simulations of IBrada fc Milgro"ml 
(|2000a[ ) show that disrupting satellites form stellar 
streams (see their Fig. 5). Thus the MOND hypothesis 

7 We do not, at present, have sufficient information to compute 
gext for the spirals in Figfl] 



for explaining deviations from the BTFR becomes very 
similar to the stellar stripping hypothesis (§ 13 . 3|) . An 
important difference is that we need not invoke highly 
eccentric orbits for the dwarfs. Since r ti M ~ r ttP hot, stel- 
lar stripping may be ongoing at present. 

The seve rity of the effect of the e xternal field can be 
quantified (|Brada fc Milgroml |20"00a) by 



7 = 



D\ 3/2 / m\!/2 



M 



(6) 



This measures how adiabatic the effects of the external 
field are. In effect, 7 is the number of orbits a star 
should make within a dwarf for every orbit the dwarf 
makes about its host. For specificity, we compute 7 for 
a star at the deprojected 3D half light radius r 1 / 2 with 
m = Mb/ '2. Equation © applies when the field of the 
host dominates: m/r 2 < Al/D 2 . Based on their bary- 
onic masses, all dwarfs are in this regim^l except Leo T. 
For th is galaxy, 7 = (D/r)(m/Af) 1 / 4 (|Brada fc MilgromI 
l2000aD . 

In order for a dwarf to remain in equilibrium, there 
must be enough internal orbits for the dwarf to adjust to 
variations in the external potential. This raises the ques- 
tion of how many orbits suffice to maintain an adiabatic 
variation. Certainly there should be at least two inter- 
nal orbits per significant change in the potential. The 
potential of the host is due entirely to the baryonic mass 
as there is no dark matter halo in MOND, so a dwarf on 
a polar circular orbit sees the host vary from face-on to 
edge-on in one quarter of an orbit. This is a significant 
variation in the potential that will only be magnified for 
non-circular orbits. Combining these criteria, we sug- 
gest that a minimum value for the adiabatic condition to 
hold is 7 > 8. This crude argument ap pears to be consis- 
tent wi th the numerical simulations of B rada fc MilgromI 
(|2000aL see their Fig. 4): satellites that remain distant 
enough from their host to always have 7 > 8 experience 
only small perturbations. Satellites reaching 7 < 8 near 
the pericenters of their orbits puff up in size and decline 
in velocity dispersion, but recover as they recede. Satel- 
lites that enter the regime 7 < 2 experience unchecked 
growth in size, subjecting them to tidal disruption and 
the formation of tidal streams. The velocity dispersions 
of these systems initially decline but then grow. This is 
qualitatively the right vector for an evolutionary track 
that would explain the deviations from the BTFR, but 
considerably more work would be required to predict 
such a track or even to check how sensitive the evolu- 
tion is to parameters such as orbital eccentricity. 

Fig. [7] shows the deviations from the BTFR and the el- 
lipticity of the dwarfs plotted as a function of 7 computed 
at the half light radius. Correlations with both are appar- 
ent. The deviation sets in around 7 ~ 10, consistent with 
the above reasoning. The most distorted dwarfs are those 
with the least time to react to changes in the external 
field. Dwarfs with large 7 have many orbits to adjust to 

8 If instead we use of os /7*i/2 to estimate the internal accelera- 
tion, all dwarfs are near the boundary gi n t ~ g ex t between domi- 
nation by the internal and external field. However, we expect the 
velocity dispersion to be inflated by the external field effect, so 
the apparent baryonic m ass is the better indicator of regime (see 
IBrada fc ATi lgrom 2000aD. 
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changes in the potential; these objects are mostly round. 
While it may be possible to avoid distorting the shapes of 
dwarfs via ti des when they are p rotected by a cocoon of 
dark matter (jMunoz et al.ll2008f). tidal distorti on seems 
inevitable in MOND (|Brada fc Milgromll2000aft . 

We note that the plot of the deviations in Fig. [7J is 
very similar to Fig. [B] One interesting difference is the 
case of Leo T, which is in the pure MOND regime. Con- 
sequently, it executes fewer internal orbits per external 
orbit than it would if it were in the quasi-Newtonian 
regime like the rest of the dwarfs. This moves it closer to 
the clump of dwarfs that adhere to the BTFR in Fig. [7J 
than in Fig.[S] We also note that Leo T is the only dwarf 
considered here that retains gas. Tides can strip gas, 
with the consequences discussed in £13.2.31 It may there- 
fore be significant that all of the remaining dwarfs are 
dominated by the external field of the host. The absence 
of gas in dwarfs within ~ 250 kpc might be a signature 
of the transition to domination by the external field in 
MOND. 

The structure of Figs. [6] and Fig. [7j should be similar 
in MOND, but this behavior is not guaranteed. It only 
follows if the MOND formula for 7 (equation |6|) provides 
the correct timescale. This would not follow had we used 
a Newtonian estimator for the timescale, either with dark 
matter or without. With dark matter, orbits within the 
dwarfs are more rapid thanks to the extra mass, and they 
always have time to equilibrate (7 ^> 10). Without dark 
matter in the dwarfs themselves, the dark matter ha- 
los of the hosts dominate and the dwarfs should dissolve 
rapidly (7 -C 1). Only when we use the baryonic mass 
for both dwarf and host in equation ^ does the onset of 
the deviation from the BTFR correspond to when inter- 
nal and external orbital timescales are comparable. This 
is one sense in which the external field effect is unique to 
MOND, and appears to be consistent with the data for 
the MW dwarf spheroidals. 

A separate question is whether the number of dwarfs 
currently near total dissolution is reasonable. By rea- 
sonable, we mean a dissolution timescale that is shorter 
than a Hubble time so that there is time for deviation 
from the BTFR to have occurred, but not so short that 
we must be catching them at a special time. Unfor- 
tunatel y, it is not clear what the dis solution timescale 
is (see ISanchez-Salcedo fc Loral I2010D . The trends in 
Fig. 5 are suggestive of progressive mass loss, but we 
do not know what the mass loss rate is, so cannot quan- 
tify t dissolve = {Mb/Mbl. As a proxy, we can estimate 
the crossing time t cross — 2r ttP hot/ ^los- Presumably it 
takes a number of crossing times for a dwarf to dis- 
solve, though it is unclear how many. Most dwarfs have 
10 8 < t cross < 10 9 years. The MW satellite with the 
longest crossing time is Sextans (~ 1 Gyr), which is a 
fair fraction of its orbital period (~ 3 Gyr if a = D). 
This object provides an interesting constraint, as it has 
had time to complete four circular orbits and a dozen 
crossing times in a Hubble time. This seems like plenty 
of time for a low 7 system so close to its tidal radius 
to have suffered some damage, and yet it remains con- 
sistent with the BTFR. If instead it is in a somewhat 
eccentric orbit that only now places it in the regime of 
tidal susceptibility, its crossing time is long enough that 
perhaps we do not (yet) see deviation in this case simply 



because the system has not had time to react. On the 
opposite extreme, the smallest of the ultrafaint dwarfs, 
Segue 1 and Willman 1, have the shortest crossing times: 
tcross ~ 5 x 10 7 years. This is uncomfortably short. Both 
are far from the BTFR, so perhaps have been caught just 
before total dissolution. Nevertheless, these cases would 
appear to pose a dissolution timescale problem. It is hard 
to say how severe this problem is without knowledge of 
the initial population. The mass function is not known 
in MOND, nor is it even cl ear whether the dw arfs are 
primordial or tidal in origin (jGentile et al.|[2007l ). 

The properties of the deviant dwarfs provide a strong 
test of the MOND hypothesis. If they are stable, undis- 
turbed systems, then their mass-to-light ratios are unac- 
ceptably large and MOND fails. If instead they are being 
tidally stripped, this situation is naturally understood in 
the context of MOND for which the tidal radii are com- 
parable to the luminous size. This differs from the stellar 
stripping hypothesis in the dark matter context in that 
tidal disruption should be going on now for the deviant 
dwarfs, and not just when the dwarfs are near the peri- 
centers of their orbits. In addition, the fact that the tidal 
radii of the dwarfs seem to be adequately estimated by 
assuming a ~ D implies that their orbits are not highly 
eccentric. The distribution of orbital eccentricities of the 
dwarfs may therefore provide another test to distinguish 
MOND from ACDM. 

4. CONCLUSIONS 

We have examined the adherence of the Local Group 
dwarf satellites of the Milky Way and M31 to the Bary- 
onic Tully-Fisher Relation. We find that most of the 
brighter dwarfs are largely consistent with the extrapo- 
lation of the BTFR fitted to isolated, late type, gas rich, 
rotating disk galaxies. The fainter dwarfs, especially the 
ultrafaint dwarfs, are not. More importantly, we find 
that residuals from the BTFR are not random, corre- 
lating well with luminosity and ellipticity. The amount 
of deviation from the BTFR also correlates with metal- 
licity, size, and the susceptibility of the dwarfs to tidal 
perturbation. We have considered a number of possible 
interpretations for the observed behavior, as we summa- 
rize below. 

Insufficient Kinematic Accuracy: — Heroic efforts have 
been made to find new dwarfs and to measure their ve- 
locity dispersions. This is a challenging endeavor. While 
the deviations of some dwarfs from the BTFR are for- 
mally significant, that significance is not overwhelming 
(typically 2 to 4a). Moreover, the analysis assumes that 
the dwarfs are spherical and in stable equilibrium. The 
assumption of sphericity at least is violated for the most 
deviant dwarfs. We therefore consider one possibility to 
be that there are no genuine deviations from the BTFR. 
This hypothesis predicts that as the data improve, so too 
will agreement with the BTFR. 

Gas Removal: — Dwarfs that deviate from the BTFR do 
so in the sense that they seem to be lacking luminos- 
ity for their velocity dispersion. This may be explained 
if baryons are removed before they form stars. Several 
possible mechanisms to accomplish this removal include 
the suppression of star formation by cosmic reionization, 
ionization from Pop. Ill stars, removal of cold gas by 
supernova feedback, and ram pressure stripping. 
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Figure 7. The residuals from the BTFR (left) and the ellipticities of dwarfs (right; symbols as per Fig. [TJ correlate with the number 
of orbits a star should complete within a dwarf for every orbit the dwarf completes about its host (equation [Sj . We expect pronounced 
non-adiabatic effects when 7 < 8 (see text). This corresponds approximately to where the discrepancy from the BTFR becomes significant, 
and to where dwarfs tend to become non-spherical. 



Cosmic reionization is often invoked in the context of 
the dwarf spheroidals, and is an attractive solution if 
only these objects are considered. If we simultaneously 
consider slightly larger gas dominated disk galaxies, it 
becomes clear that reionization is not in itself an ade- 
quate explanation for the observed trends in the data. 
Some other mechanism must be acting to suppress the 
accumulation of cold baryons in a manner that becomes 
more severe with decreasing V c . Whatever this mecha- 
nism is, it presumably affects the dwarf spheroidals as 
well. Cosmic reionization may be an additional factor 
acting only at scales V c < 20 kms -1 . 

Feedback from supernovae is a candidate mechanism 
for affecting star formation across all halo masses. This 
provides a qualitatively appealing explanation for the 
trend in the detected baryon fraction with halo mass. 
Supernovae provide the kinetic energy to drive gas be- 
yond escape velocity, but the escape velocity increases 
with increasing halo mass so progressively more baryons 
are retained. This mechanism may provide a natural ex- 
planation for the residual correlation with luminosity and 
metallicity. However, quantitative tests remain wanting, 
as does an explanation for the correlations of the residu- 
als with ellipticity and tidal susceptibility. 

In the ram pressure stripping hypothesis, galaxies de- 
viate from the BTFR when they fall into the halo of 
the current host galaxy and their cold gas is ablated by 
the ram pressure of the hot gas in the halo. This hy- 
pothesis requires that a sufficient amount of hot gas be 
present in the halos of the host galaxies, which is not ob- 
viously the case. It predicts that the ages of stars in the 
dwarfs is related to the time of infall, with the dwarfs 
that fall in first losing the most gas, deviating by the 
largest amount from the BTFR, and having the oldest 
stars. Strictly speaking, the star formation history of 
the pre-infall dwarf is not constrained, so there should 
be some scatter in these predictions. However, no star 
formation can occur after infall and gas stripping, so the 
ages of the youngest stars should correspond well to the 
time of infall with a sharp truncation in star formation 
after that time. This would provide a natural explana- 
tion for why some dwarfs appear to have had relatively 
recent star forming events but now contain no cold gas. 



One would also expect the metallicities of the stars to 
reflect the star formation history. The first dwarfs to 
fall in would have had the least time for enrichment and 
have the lowest [Fe/H]. This predicts that [a/Fe] should 
also correlate with the amplitude of deviation Ft,, to the 
extent to which we expect the objects with the briefest 
enrichment time to have the highest [a/Fe]. 

The gas removal hypotheses provide a potential expla- 
nation for the correlation of BTFR residuals with lumi- 
nosity and metallicity. However, they provide no obvious 
explanation for the correlation with ellipticity and tidal 
susceptibility. This occurs more naturally in the follow- 
ing two scenarios. The removal of gas and the concomi- 
tant truncation of star formation and its consequences 
for metal enrichment may also occur as a result of tidal 
stripping in the following hypotheses. 

Stellar Stripping: — The correlation of the residuals with 
ellipticity in addition to luminosity suggests a role for 
tidal effects. In this hypothesis, the dwarfs become pro- 
gressively more distorted due to tidal disruption as they 
orbit the Milky Way. Stars are lost in the process, reduc- 
ing the luminosity of the dwarfs. At the present time, the 
computed tidal radii of the dwarfs greatly exceed their 
luminous extent. This leads us to infer that the orbits of 
the dwarfs must be highly eccentric in this scenario, with 
most of the stripping occurring during pericenter pas- 
sage. This hypothesis predicts that the mass required to 
reconcile each dwarf with the BTFR may exist in a tidal 
stream. It further predicts that the age and metallicities 
of stars in the predicted streams should be consistent 
with those of the parent body. Examples exist where 
this may already be observed. 

MOND: — Low surface brightness dwarf spheroidals pro- 
vide a strong test of an alternative to dark matter, 
MOND. They should very nearly follow the BTFR, which 
is a consequence of the specific form of the modified force 
law in MOND. While some dwarfs do indeed adhere to 
the BTFR, others deviate substantially. The ultrafaint 
dwarfs of the Milky Way have MONDian mass-to-light 
ratios in the tens to hundreds. This is fatal for MOND 
if these dwarfs are in a stable equilibrium, the spheri- 
cal approximation used in the analysis is adequate, and 
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the kinematic data are to be trusted. Intriguingly, the 
sizes of the dwarfs relative to their MONDian tidal radii 
correlate strongly with the degree of deviation from the 
BTFR. Indeed, the discrepancy for MOND sets in pre- 
cisely where the theory predicts that non-equilibrium ef- 
fects become strong. It therefore appears that the unac- 
ceptably high mass-to-light ratios may be a result of the 
dwarfs being out of equilibrium. This should be testable, 
in the sense that the deviant dwarfs should show evi- 
dence of tidal disruption while the dwarfs that adhere to 
the BTFR should not. Notably, stripping of the deviant 
dwarfs should be ongoing and not restricted to pericenter 
passage as in the stellar stripping hypothesis. 

It is of course possible that some combination of these 
effects is at work. As a dwarf satellite approaches its host 
on its orbit, it is subject to both tidal forces and ram 
pressure effects. Perhaps gas is lost first due one or both 
of these effects, with stars being tidally liberated later. 
This makes it somewhat difficult to distinguish between 
the various hypotheses, but it should be possible. 
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